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Abstract: Input from QCD and string theory is used in order to elucidate basic 
features of the string theory dual of QCD, It is argued that the relevant string 
theory is a five- dimensional version of the type-0 superstring. The vacuum solution 
is asymptotically AdSs, and the geometry near the boundary is stringy. The 
structure of YM perturbation theory however emerges near the boundary. In the IR, 
the theory is argued to be well-approximated by a two-derivative truncation that 
takes into account strong coupling effects. This explains the success of previously 
proposed five- dimensional Eistein-dilaton gravity with an appropriate potential to 
describe salient features of the strong YM dynamics. 
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1. Introduction and outlook 

QCD is a very successful theory of strong interactions. It is also a theory that is 
hard to calculate with, due to the strong coupling region in the IR. As any kind of 
observable physics passes via low-energy filters it has complicated efforts in the past 
three decades to test the theory and make predictions. Our theoretical understanding 
of QCD stems from several sources/techniques. 
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• Weak-coupling perturbation theory. This is valid in the UV, because of asymp- 
totic freedom and has been a key element in the understanding of the strong 
force. Its impact is relying of the factorization of processes into hard and soft 
components. 

• Euclidean Lattice techniques based on numerical estimates of the QCD path in- 
tegral. To date this is the most direct non-perturbative technique that has pro- 
vided the first clues to confinement, and numbers for the hadron spectrum that 
could be compared with data, among other things. This is a non-perturbative 
approach that is fully mature and its main technical limitation today is com- 
puting power. As it is inherently Euclidean it cannot however address ab initio 
a class of problems that describe time-dependent dynamics. This includes fi- 
nite temperature dynamical phenomena, as well as scattering. Although some 
quantities can be obtained by analytic continuation extra input is needed in 
order for analytic continuation to be performed reliably. 

• Special purpose phenomenological models and approximations. For specific 
problems, phenomenological models can give deep insights into physics that is 
not directly under analytical control in QCD. One notable example is Chiral 
Perturbation Theory. This is a low-energy effective field theory for the light me- 
son sector based on ideas of chiral symmetry breaking. Other examples include 
resumations of perturbative effects based on some assumptions, truncation and 
solutions of Schwinger-Dyson equations, applications of the Large- N c expansion 
and associated matrix models, etc. 

In the past decade there have been two developments that stirred the field of 
strong-interaction physics. The first is data from the RHIC collider that gave the 
first solid indications for the physics of the quark-gluon plasma Q • The second is new 
intuition and results on the large- N c expansion of gauge theories that changed our 
perception of the description of strongly-coupled large- N c gauge theories, 0, 0, |]. 
The prototype example has been the AdS/CFT correspondence as exemplified by the 
(well studied by now) duality of Af = 4 super Yang-Mills theory and IIB string theory 
on AdS$ x S 5 . Further studies focused on providing examples that are closer to real 
world QCD, j5], [|. It is fair to say that we now have a good holographic understanding 
of phenomena like confinement, chiral symmetry and its breaking as well as several 
related issues. The finite temperature dynamics of gauge theories, has a natural 
holographic counterpart in the thermodynamics of black-holes on the gravity side, 
and the thermal properties of various holographic constructions have been widely 
studied, 0, [7], H || [H|, exhibiting the holographic version of deconfinement and 
chiral restoration transitions. 

The simplest top-down string theory model of QCD involves D 4 branes with 
supersymmetry breaking boundary conditions for fermions 0], as well as a flavor 
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sector that involves pairs of D$ — D$ probe branes inserted in the bulk, [II]]. The 
qualitative thermal properties of this model closely mimic what we expect in QCD, 
[§]. Although such theories reproduced the qualitative features of IR QCD dynamics, 
they contain Kaluza-Klein modes, not expected in QCD, with KK masses of the same 
order as the dynamical scale of the gauge theory. Above this scale the theories deviate 
from QCD. Therefore, although the qualitative features of the relevant phenomena 
are correct, a quantitative matching to real QCD is difficult. 

Despite the hostile environment of non-critical theory, several attempts have been 
made to understand holographic physics in lower dimensions in order to avoid the 
KK contamination, based on two- derivative gravitational actions, Jl2, 13]. Indeed, 
large N QCD is expected to be described by a 5-dimensional theory. The alternative 
problem in non-critical theories is that curvatures are of string scale size and the 
truncation of the theory to the zero mode sector is subtle and may be misleading. 

A different and more phenomenological bottom-up approach was developed and 
is now known as AdS/QCD. The original idea described in [TJ[] was successfully 



applied to the meson sector in fLo| , and its thermodynamics was analyzed in [gj. 
The bulk gravitational background consists of a slice of AdSs, and a constant dilaton. 
There is a UV and an IR cutoff. The confining IR physics is imposed by boundary 
conditions at the IR boundary. This approach, although crude, has been partly 
successful in studying meson physics, despite the fact that the dynamics driving 
chiral symmetry breaking must be imposed by hand via IR boundary conditions. 
Its shortcomings however include a glueball spectrum that does not fit very well the 
lattice data, the fact that magnetic quarks are confined instead of screened, and 
asymptotic Regge trajectories for glueballs and mesons that are quadratic instead of 
linear. 

A phenomenological fix of the last problem was suggested by introducing a soft 
IR wall, Although this fixes the asymptotic spectrum of mesons and meson 

dynamics is in principle self-consistent, it does not allow a consistent treatment of 
the glue sector both at zero and finite temperature. In particular, neither dilaton nor 
metric equations of motion are solved. Therefore the "on-shell" action is not really 
on-shell. The entropy computed from the BH horizon does not match the entropy 
calculated using standard thermodynamics from the free energy computed from the 
action, etc. 

A well-motivated way to obtain linear Regge trajectories for mesons was advo- 



cated in In particular it was pointed out that the natural order parameter for 
chiral symmetry breaking in the context of flavor branes is the open-string tachyon. 
By studying the tachyon dynamics, it was shown that confinement in a wide class 
of backgrounds is enough to guarantee chiral symmetry breaking, with linear meson 
trajectories. When flavor branes are at distances larger than the string scale, the 



analogue of the tachyon order parameter was investigated in ||22|| . This is relevant 
for flavor sectors that resemble more the Sakai-Sugimoto setup. 
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An improved holographic model that lies somewhere between bottom up and 
top down approaches has been proposed, [17], [18| . It is a five-dimensional Einstein 



dilaton system, with an appropriately chosen dilaton potential. The vacuum solution 
involves an asymptotically logarithmically AdS solution near the boundary. The bulk 
field A, dual to the 't Hooft coupling, is vanishing logarithmically near the boundary 
in order to match the expected QCD behavior. This implies that the potential must 
have a regular Taylor expansion as A — > 0, and that A = is not an extremum of 
the potential. This is unlike almost all asymptotically AdS solutions discussed so 
far in the literature. In particular the canonically normalized scalar (the dilaton) 
is diverging at the boundary r — > as <f) ~ — log(— logr). The coefficients of the 
UV Taylor expansion of the potential are in one-to-one correspondence with the 
holographic /^-function. 

In the IR, the potential must have an appropriate behavior so that the theory 
is confined, has a mass gap and a discrete spectrum. This selects a narrow range of 
asymptotics that roughly obey 

V(X) ~ \ 2Q , A -> oo with ~ < Q < -. (1.1) 

o o 

The vacuum solution always ends in a naked singularity in the bulk. Demanding 
that this is a "good" singularity in the classification of Gubser |23| implies Q < 
4/3. Simple interpolations between the UV and IR asymptotics reproduce very well 
the low-lying glueball spectrum as well as the perturbative running of the 't Hooft 
coupling fI5|| . At finite temperature this model exhibits the behavior expected from 



QCD: There is a deconfining transition and the thermodynamics is very close to what 
one expects from lattice QCD, |05], ^(J. 



In this paper we will go through several arguments originating both in string 
theory and QCD as we understand it, that will help us analyse in more detail the 
structure of the string theory dual of QCD. We will see that although there are 
ambiguities in several places, a picture emerges that seems consistent and gives some 
hope that we may one day tame the associated string theory. Even today, it may be 
used as qualitative litmus test of ambitious holographic models. 

There are several directions that have not yet been explored. An important one 
concerns the behavior of one and two-point functions in the IR. This is an important 
issue as it stands at the heart of justifying the neglect of vevs of higher-dimension 
operators in a holographic context. Although techniques similar to what we use can 
be used in this direction, we will not attempt this here. 

Another issue is the theoretical definition and practical viability of a hybrid 
model for QCD. In such a model, physics in the UV is described via perturbative 
QCD that is used to generate boundary conditions, at a rather low scale (in the 
few GeV region). Below this scale a holographic model should be used. IN such a 
hybrid model, the UV region near the boundary (that as we argue here is stringy) 
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can be altogether avoided. The IR region (that as we argue here can be reasonably 
well-described by a two derivative action) can be handled with standard holographic 



techniques. An attempt in this direction can be found in . 

An interesting issue is the cosmological evolution of strongly coupled matter, 
both made of glue and quarks. In the former case the setup is almost identical to 
the one studied in the context of Randall-Sundrum cosmology as was shown in P5|] . 
Indeed the simple solution for conformal matter described in has an alterna- 
tive description in terms of lowering the UV cutoff brane inside AdS. The Randall- 
Sundrum tuning corresponds to the choice of coupling gravity to renormalized rather 
than bare sYM operators (the vacuum energy in particular). A similar study for a 
non-conformal theory like QCD involves a few extra ingredients, the most important 
of which is establishing the geodesic motion of the UV boundary in the bulk, and in 
particular the dilaton couplings to the boundary. This is interesting as it may give 
new tools to study the impact of the deconfinement phase transition in a cosmological 
setup. 



2. General remarks on the string theory dual 

The first question we may pose is: in how many dimensions is the string theory dual 
of QCD living? A way to answer this question proceeds via the intuition developed 
in the past 20 years from matrix and other large N theory duals to string theory. 
Indeed, the intuition is as follows. The large N-gauge theory contains several adjoint 
fields living on a d-dimensional space M^. Typically the eigenvalues of the adjoint 
matrices becomes new continuous dimensions. For example in the case of the "old 
matrix models" the single eigenvalue distribution increases the spacetime dimension 
by one. 

Not all adjoint fields provide independent eigenvalue distributions and therefore 
holographic dimensions. Fields related by global symmetries must be reduced appro- 
priately. Also if the symmetry is not there, but there is a symmetric theory related 
by RG flow to the previous one, than there is a reduction in the number of holo- 
graphic dimensions. For example in M = 4 superYM, there are 4 adjoint matrices 
from the vectors, 8 adjoint matrices from the fermions and 6 adjoint matrices from 
the scalars. Four-dimensional Lorentz invariance of the vacuum indicates there is 
a single independent eigenvalue distribution from the vectors. Similarly the SO (6) 
symmetry implies that that there are 5 independent eigenvalue distributions from 
the scalars (nicely exposed in the work of Berenstein and collaborators, |26|| ) 

In the case of QCD the situation is simpler. We have four (a vector) adjoint 
fields. Since the theory and its vacuum on flat space are expected to be Lorentz 
invariant, only one eigenvalue distribution is independent. Therefore we expect one 
extra holographic dimension and therefore the string theory dual is expected to live 
in 5 dimensions. 



- 5 - 



QCD unlike successful gauge theories in the holographic domain, is asymptot- 
ically free. This is a property that has marred attempts in the past 10 years to 
work out a trustworthy dual string theory. Intuition coming from M = 4 superYM 
indicates that the 't Hooft coupling is directly related to the spacetime curvature as 
A ~ -^p-. If this is taken at face value, it would imply that the putative QCD dual 
will have singular curvatures near the UV, where the QCD coupling vanishes. 

There are several caveats to this line of reasoning, that we will mention here. 

• The relation between the 't Hooft coupling and spacetime curvature as implied 
by M = 4 superYM may not be universally applicable and in particular its ex- 
trapolation to weak coupling may not be trustworthy. We have many examples 
of such behavior, where non-linear behavior of effective actions (notably the 
DBI and CFT actions) smooths out the leading singular behavior. Controlable 
examples in closed string theory itself, that include WZW and coset models do 
indeed behave differently at strong curvature. A suggestive example is SU(2)fc 
WZW model where the limit of strong curvature has either curvature of order 
the string scale (k=l) or contains no space at all (k=0). 

It is plausible that a respectable alternative is that in the limit of the vanishing 
coupling, curvatures remain at the string scale. Therefore although the regime 
is stringy, it is not singular. 

• The 't Hooft coupling in J\f = 4 is constant. As such, one can rescale the metric 
and compute the curvature in different "frames". Although different rescaled 
curvatures behave differently, going from one to the other is a simple process. 
On the other hand the choice of frame is a relevant question in the string 
theory dual of QCD, where A is expected to be a function of the holographic 
coordinate, reproducing the RG running of the YM coupling. Different frames 
with include derivatives of the A in the curvature and can radically alter its 
behavior. 

• Another issue is the approximate conformal invariance that characterizes QCD 
in the extreme UV. We do expect that this will be geometrically encoded in the 
string theory dual. The generic lesson of holography in controlled examples is 
that conformal invariance and its approximations are encoded in the structure 
of the string geometry via AdS spaces and their deformations. AdS is a very 
symmetric spacetime, and carries a single boundary where the UV definition 
of the holographic theory resides via sources. It would be a rather unfortunate 
situation that the space is singular at the place where we are supposed to define 
the UV of the theory. 

Because of all of the above, a reasonable expectation of the structure of the string 
background describing the QCD dual is a metric that is asymptotically AdS$ near its 
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boundary region (UV). We expect that this will be modified as we flow towards the 
IR signaling the breaking of conformal invariance. Moreover as we will argue more 
in detail further on, the AdS 5 curvature near the boundary should be comparable to 
the fundamental string length £ s . 

2.1 The low-energy string spectrum: a gauge theory view 

An important question that is crucial in setting up the vacuum problem in the string 
theory dual is to estimate which of the bulk fields are important in determining the 
vacuum structure. The lowest dimension fields are certainly the most important in 
the UV. QCD has no (strongly) relevant operators. The first non-trivial operators 
start at dimension 4 and are given by the quadratic trace of the field strengths 



Tr[F^F pa ] 

To distinguish different operators we may use the U(d) decomposition 



(B®B 



symmetric 



(2.1) 



(2.2) 



Finally we must remove traces to construct the irreducible representations of 0(d): 



>me« , 



(2.3) 



where the line under a Young tableau implies that all possible traces have been 
removed and it therefore represents an irreducible representation of 0(d). A • stands 
for the singlet of 0(d). 

The two singlets are the scalar (YM Lagrangian, dual to the string theory dilaton) 
and pseudoscalar (instanton) densities: 

(j)^Tr[F 2 ] , a <-> Tr[F A F] (2.4) 

Each carries a single d.o.f. The t' Hooft coupling A is related to the dilaton as 

A ~ (2.5) 

at least in the UV. The instanton density should be dual to a pseudoscalar bulk 
axion, a. The next operator is a traceless conserved symmetric tensor 



m 



Tpv = Tr 



, (T)/ = , d^ u + 



(2.6) 



It is dual to the 5-dimensional graviton, g^ v . Finally 



^ iiv\po ~ ^ r \F nvF pu ^{d^pFy^ QvpF^ 9fiaF L , p +g ua F^ p ) + ^ {g^p9va 9vp9po)F 

(2.7) 



-7- 



It has 10 independent d.o.f and should be dual to a similar massive tensor in the 
bulk string theory. 

Unlike M = 4 superYM, near the UV, the dimensions of these operators are 
reliable, as they are given by their free-field theory values plus small corrections. Of 
all these operators with A = 4, only the last one correspond to a field that is massive 
in the bulk string theory. Therefore it is expected to be less important at least in 
the UV. 

Operators with A = 5 are given by Tr[V M F vp F (TT ] (with V being the gauge 
covariant derivative) while operators with A = 6 are given by Tr[V ^F^pV ^ F aT ] and 
TrfF^F^vFpo-]. Of all the higher dimension operators one is expected to correspond 
to the NS antisymmetric tensor, namely 

~ Tr[F [)ia F ah F bu] + l -F ab F ah F^] (2.8) 



[34]] and should be massive with a UV mass M = j^— in order to have the correct 
scaling dimension in the UV. This is happening because of a combination of two 
effects : the fact that B appears in the five-form field strength together with the RR 
two-form C*2 and the fact that the RR-five form has a vev. This provides a mixing 
of B 2 and C 2 that gives both a mass. This is similar to what happens in Af = 4 
superYM. 

The conclusion of this section is that the massless bulk fields should be dual to 
the operators, Tr[F 2 ] (dilaton), T^ v (metric), Tr[F A F] (axion). 

2.2 Bosonic string or superstring? 

We will now present an argument that suggests that the string theory dual to QCD 
should have a RR sector, which furthermore implies that it is a superstring theory. 
This does not necessarily imply (broken) spacetime supersymmetry but rather that 
the world-sheet gauge symmetry of the theory is some form of supersymmetry, which 
implies in particular the existence of a (bi-fermionic) , RR sector. Note however that 
pure YM does not have gauge invariant fermionic operators. Therefore, the string 
theory, although a superstring theory, it should not contain spacetime fermions (NS- 
R, and R-NS sectors). Type theories, have precisely this property, and have been 
candidate dual grounds for YM for some time |28[ . 



It should be noted here that even after we add quarks, all gauge invariant opera- 
tors in QCD (with the exception of baryon operators) are bosonic. Baryon operators 
on the other hand should correspond to appropriate "solitonic" D-branes in the string 



theory as we now understand in many similar examples |3(J. Therefore the standard 
spectrum of the string theory dual to QCD should also contain no spacetime fermions. 

A first candidate for a RR field should be the RR four-form, C4, that in standard 
examples is known to provide the flux responsible for introducing the (large number 
of ) -D3 branes into the background geometry. Five dimensions is a very special 
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dimension however for C4 as it does not contain any propagating degrees of freedom, 
and this makes therefore its presence a bit murkier. There is a however another 
bulk field that (a) is propagating and (b) sources one of the couplings of the YM 
theory. This is the axion, dual to the YM instanton density. This field must be a RR 
field (as indeed happens in the TV = 4 superYM example), in order to match known 
properties of the CP-odd sector of large- N c YM []27J • 

The action of large- N c QCD including a 9 angle can be written in the form 



S- 



YM 



d x Tr 



—F 2 + 

4g* 



9 



-F A F 



N r / <Tx Tr 



F 2 + 



9 



■F A F 



(2.9) 

The proper scaling in the 't Hooft large-iV c limit is to keep A = g 2 N c and ( = jj~ 
finite and fixed. Consider now the ^-dependent vacuum energy, to leading order in 
the 1/N C expansion 

E YM = N 2 F[\X} (2.10) 

This must be invariant under the #-angle periodicity shift 9 — > 9+2-k. This however is 
impossible if F is a smooth function. This can be achieved only if F is multibranched 
function (obtained by minimizing a collection of many nearly degenerate minima). 
In view of this we can write the vacuum energy as f27|| , 



E YM = N 2 c Min keZ f 



A. 



9 + 2nk 



(2.11) 



which shows that it is periodic but it is not a continuous function of 9. The CP 
transformation 9 — > —9 implies that f(X,9) = f(X,—9). CP is unbroken only if 
9 = 0,71. 

The integer k labels different vacua that are related by integer shifts of the 9- 
angle. The absolute minimum is expected at 9 = at the k = vacuum. Taking all 
this into account we can write at large N c 



E YM = N 2 f[X,0] + 



1W] 
2 d( 



Min kez (9 + 2ixk) 2 + 0(N~ 2 ) 



(2.12) 



The quantity x = §^ [A, 0] is known as the topological susceptibility and it is known 
to be non-zero |31]. We therefore observe that the leading 9 dependence of the 
vacuum energy is coming in at order 0(1) while higher terms are further suppressed 
with N c . 

We will now show that this property on the string theory side is due to the 
special properties of RR fields. The axion is dual to the instanton density. This 
implies that its source gives the UV value of the 9 angle 



a(r) = 9 UV + 0{r 4 



(2.13) 
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Assuming it is a RR field we may write its tree-level effective action as 



S = M 3 J (fxyfge 



2,p 



R + 4(<90) 2 + ^-(da) 2 + C A e^(da) i + 



(2.14) 



where it should be noted the peculiar dependence on the dilaton of the axion terms 
in agreement with standard string-theory dilaton counting. We now translate to 
variables that have a smooth large-iV c limit A = N c e^ to rewrite the action with 
explicit N c dependence 

2 



S = N 2 M 3 j d 5 x^g± 



(2.15) 



As in the on-shell action (da) 2 ~ 6 2 , (da) 4 ~ 8 A etc, we observe that we obtain the 
same large- N c scaling of the different ^-dependent terms as in the field theory side , 
( |2.12| ). More details on the vacuum solutions and action for the QCD axion can be 
found in |T§. 

The upshot of the previous analysis is that the axion in the string theory dual 
of QCD is a RR field, indicating that the string theory is a superstring theory in the 
type-0 class. 

2.3 The minimal low-energy spectrum: a string-theory view 

From our discussion so far we have argued that the string theory dual must have a 
NS-NS sector with the usual fields, g^, B^, 0, as well as a RR sector that contains at 
least the axion a = C and the four-form gauge potential C 4 necessary for generating 
the color flux, which will be proportional to the number of colors iV c . 

Two main issues require a discussion. The first is the RR sector. The minimal 
possibility is a spinorx spinor in five dimensions. Five-dimensional spinors are not 
chiral, therefore no gaps are expected in the expansion of the RR bispinor. A direct 
group-theoretic expansion gives 



spinor 5 x spinor 5 = F + F 1 + F 2 + F 3 + F 4 + F 5 



(2.16) 



where F p = dC„_i is a p-form field strength. Moreover, the truly independent fields, 



are ^0,1,2 as the rest are related by Poincare duality: F 5 
as implied by the properties of spinors. 



"Fn, Fa 



*Fi, F 3 



F 5 = *F . F 5 ~ N c and the associated C4 is not a propagating field in five 
dimensions. Therefore the 5-flux is important as a background, and as argued 



in []17 , p~8|] is responsible for non-trivial dilaton dependence in the tree-level 
string effective action. It provides an IR effective potential for the dilaton 
(=QCD coupling), and other important dynamical effects in the flavor sector, 
in particular the correct flavor anomaly generating CS terms [El], H- Having 



no propagating degree of freedom, it is not dual to any operator of large- N c 
YM. 
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i*4 = *F\. F\ is the field strength of the axion a dual to the QCD instanton 
density. Its dual form is a three- form C 3 . 

F 3 = *F 2 . F 3 is the field strength of a RR two-form, G 2 . Its dual is a one- 
form C\. Although these forms belong to the lower level of the RR spectrum, 
they should correspond to massive states of the string theory, and therefore to 
higher (than four) dimension operators in QCD. C 2 becomes massive due to 
its (CS related) mixing with E>2 from the NS-NS sector. Ignoring the axion we 
may write the leading order (string-frame) action for the three- and five-forms 
as 



S = — M 3 J d 5 x^/g 



(2.17) 



2-3f 3 2-3! 3 2-5! 5 . 

F 3 = dC 2 , H 3 = dB 2 , F 5 = dC 4 -C 2 AH 3 (2.18) 

The equations of motion that stem from this action are 1 

V^(e- 2 ^ 3 ,^p) + ^upafrFffr = , WF 3 ^ p + l -F^ paPl H^ = 

(2-19) 

V M F 5 uvpar = — > F^ avpaT = M "^ T — ~n ~ (2.20) 

where w is a dimensionless constant. Substituting the five-form in the three- 
form equations we obtain, 

^(e-^H 3 ^ p ) + ^ € -^F 3 ^ = , V^ p + ^^^ = 

(2.21) 

We may decouple the two equations by direct manipulation to obtain 

W N 2 

V m [y^e- 2 ^,^ + cyclic] + ^-^H 3fUptr = (2.22) 

and a similar one for F 3 . The N c dependence is also appropriate as the equa- 
tion (|2.22j ) is iV c -independent if we introduce the appropriate ^-independent 
variable A = N c e^ proportional to the YM 't Hooft coupling in the UV. 

The upshot of this analysis is that both B 2 and C 2 combine to a massive two- 
tensor, that should be dual to the C—odd non-conserved operator Tr[F^ a F ab F bu ] + 
\F ab F ab F^} with UV dimension 6. 

Therefore, this minimal spectrum includes all fields we expect to be massless in 
five dimensions (g^, <f), a) as well as the flux-generating four-form. In simple type-0 



1 These equations are consistent only if the dilaton is constant, which is the case in the Af = 4 sYM 
case. Here non-linearities are important as we will see later one, but the basic mixing mechanism 
is similar. 
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vacua in ten dimensions, there is a doubling of RR fields. This is due to the fact 
that there is effectively no chirality projection. However, in 5 dimensions, there is 
no chirality and it is expected that the RR sector would have the form advocated 
above. 

There is another point that needs discussion: in type-0 vacua in ten and six 
dimensions there is a closed-string "tachyon" scalar. It is indeed a tachyon in flat 
ten-dimensional space [29|, but it may be massless or even massive in curved non- 
critical backgrounds [j32| . There is certainly no place for a tachyonic or massless 
scalar near the boundary of AdS in the dual string theory of YM as that would 
imply the presence of another relevant or marginal operator in YM. However all 
such operators have already accounted for. Therefore, if the string theory has such 
a zero-th level scalar, it should be massive. Moreover, it is not at all obvious what 
operator should be the dual of such a scalar 2 . 

Of the minimum set of fields that we mentioned above, only g^, <p, a, F 5 can 
have vevs (non-trvial profiles) in the vacuum in order to preserve 0(1,3) Lorentz 
invariance. In particular this precludes vectors and two-index antisymmetric tensors 
from obtaining a vev. 



2.4 The relevant charged defects 

Several strings and branes are expected to exist as solitonic objects in this string 
theory, in analogy with critical string theory. We will enumerate them trying to 
elucidate the nature of each defect. 



• The fundamental string. It couples electrically to B^ v . This is expected 
to represent the YM flux tube. Its tension, Tp ~ -m is iV c -independent. It 
should not be confused with the QCD string tension a that is multiplying the 
linear term in the inter-quark potential. This is proportional to Tp but the 
string-frame scale factor also enters fTBfl . 

• NSq brane. In five dimensions it couples magnetically to B^ u . It is the ana- 
logue of the NS$ brane of critical string theory. It is a "point-like" soliton with 
tension that scales as T^s a ~ 0(N%). It should be thought of as a magnetic 
baryon vertex that binds together N c magnetic quarks (each having a mass 
that scales as 0(N C )) . 

• D-i branes. They are the YM instantons. They couple electrically to the 
axion. Their "tension" is 0(N C ) in agreement with what we expect from the 
instanton action. 

2 Some intriguing observations on the couplings of such a scalar to probe D-branes were made in 
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• Dq branes. They couple magnetically to Ci- They are the baryon vertices in 
five dimensions. Their WZ couplings are responsible for binding N c fundamen- 
tal strings. Their tension is 0(N C ). On flavor branes they couple to baryon 
number that is equivalent to instanton number of the flavor U (Nf) gauge fields. 

• D\ branes. They are the magnetic strings, namely flux tubes between mag- 
netic quarks. They couple electrically to C 2 , and have a tension of order 0(N C ). 

• D 2 branes. They couple magnetically to the axion. They are domain walls in 
4 dimensions that separate different oblique confinement vacua. As one moves 
across a D 2 brane, 9 jumps by 2n. 

• D 3 branes. They couple electrically to the four-form. They generate the gauge 
group of the gauge theory. 

• -D4 branes. They are space filling branes that generate flavor in the YM theory. 
Since the string theory is oriented, tadpole cancelation implies that an equal 
number Nf of D4 and D4 branes must be introduced in 5 dimensions. The 
strings stretched between D 3 and D 4 branes generate the left-handed quarks 
while the ones stretched between D 3 and D A branes generate the right-handed 
quarks. 

2.5 Why a spectrum truncation might work for the vacuum structure? 

Large-iV c YM has an infinite number of single-trace operators. The UV definition 
of the theory, involves only T^ u , and Tr[F 2 }. If 9 UV is non-zero, then Tr[F A F] 
is also involved. In the holographic dual this translates into the statement that in 
the vacuum solution, only 0, g^ v and potentially a have a source term in the UV 
boundary. This guarantees that their profile in the string-theory vacuum solution 
is necessarily non-zero. Although all the other (infinite tower) of bulk fields have 
no UV sources, this does not necessarily imply that their profiles in the vacuum 
solution vanish. They can have non-trivial vevs, that would trigger a non-trivial 
solution profile in the holographic direction. 

In a theory with exact conformal invariance, non-zero one-point functions, can 
be redefined to zero by a subtraction. Because the theory is conformal, once they 
are set to zero at a given scale, they remain zero at all scales. This is not the case 
in theories where conformal symmetry is broken , as in QCD. Although one can 
subtract a one-point correlator a a given scale, this does not guarantee than the one- 
point function remains zero at all scales. Of course, bulk fields that break Lorentz 
invariance cannot acquire vevs, and therefore have trivial profiles in the vacuum. 
But this is not the case for example for bulk scalars. One therefore may ask, can 
we neglect the non-trivial profiles of all other single trace operators in the YM/QCD 
string vacuum solution? 
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The answer in the UV is rather simple: the higher the scaling dimension and spin 
of an operator, the larger its bulk mass, and the smaller its influence in the equations 
of motion of basic fields, g^ u , 0, a. This is a well-known effect both in asymptotically 
AdS space-times and asymptotically free QFTs. This is the reason we can truncate 
the infinite coupled system of RG equations near a free field theory, and study a small 
number of flows corresponding to the most relevant operators. Therefore in the UV 
of YM, (free-field) scaling dimensions and spin determine the relative importance of 
operators. 

The situation in the IR is more complicated and strictly speaking beyond con- 
trol in QCD. Generally speaking, higher- dimension operators can and do sometimes 
become important in theories that are strongly coupled in the IR. A typical example 
involves (massive) KK fields in the bulk that their vevs can be important in resolv- 
ing IR singularities in the bulk. Is this expected to happen in QCD? The direct, 
short and honest answer is : we do not know. There are however a few tantaliz- 
ing arguments and indications in the past two-three decades that point towards the 
following answer: that for many, IR relevant and simple observables, vevs of higher 
dimension (in the UV) operators are not that important for the IR physics. A large 
class of such arguments are summarized by the surprising successes of the SVZ sums 



rules |33| that seem to hint at the previous statement. A holographic argument in 
the same direction, elaborated in the context of the string theory dual will also be 
presented later on in this article. In view of this we will entertain the possibility that 
for several observables it is enough to consider the vacuum structure as described by 
"light string spectrum" , g^ u , <fi, a together with the the non-propagating four-form 
potential. 

2.6 The vacuum solution ansatz 

We will consider YM and QCD on Minkowski space, or Euclidean four- dimensional 
(flat) space. More exotic geometries maybe considered but we will not explore them 
here to keep the discussion simple. The theory is Lorentz invariant and we do not 
expect that Lorentz invariance will be broken in the quantum theory. Therefore the 
form that the metric and other fields will take in the vacuum is very constrained 

ds 2 = g ta ,dafdx v = b(r) 2 [dr 2 + dx M dx M ] , <P -> <t>(r) , a -> a(r) (2.23) 

where g^ v is the five-dimensional metric, x M are the four- dimensional Minkowski co- 
ordinates and r is the radial coordinate. One can still perform radial reparametriza- 
tions, which have been used to bring the metric to the form above. There is also the 
five-form field strength F§. 4d Lorentz invariance implies that 

^1^2— MS = f( r ) e MiM2— (2.24) 

The fact that we are in five dimensions and the fact that Lorentz-invariance con- 
straints the metric as in ( |2.23| ) already implies that capped geometries, like the one 
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in |kJ that have been popular in order to describe backgrounds that are confining in 
the IR are not possible here. The reason is that there is no extra holographic coor- 
dinate beyond the radial one to generate the cigar geometry familiar from Euclidean 
black holes. Lorentz invariance on the other hand prohibits the used of a Minkowski 
coordinate to that effect. 

From now and for the rest of this paper we will neglect the axion. As discussed in 
section 2.2 , its kinetic term is large- N c suppressed and it does not therefore contribute 
to the vacuum structure at leading order in 1/N C , except for CP-odd observables. 
The YM axion has been discussed in [O, p|. 



3. The string effective action 

In this section in preparation for our exploration of the "vacuum" of the string theory 
dual to YM we will investigate the general action at string tree level. As we will see, 
at least in the UV, the geometry is expected to be stringy, therefore we do not expect 
a few terms in the derivative expansion to be a good guide. In the absence of an 
exact string description we will use the language of arbitrary functionals of local 
curvature invariants as a guide. Although this intuition may sometimes fail we will 
trust it to obtain qualitative conclusions. 

The tree-level string effective action that we will start (in the string frame) is 

5 tree = M 3 | C / 5 x v ^e- 2 ^4(^) 2 + F( J R,0] , £ = -e 2 ^| (3.1) 

In this action we have included the fields that will be non-trivial in the vacuum 
solution, namely g^ v , 4>, F 5 . We have suppressed most possible distinct tensor struc- 
tures under which the curvature and the five-form field strength enter in the above 
action. In particular, the effective action is expected to be a non-linear function of 
the Riemann tensor and its covariant derivatives. Although such terms are impor- 
tant ingredients in several string theory observables 3 , for the arguments that will be 
made, a simplified action involving only the scalar curvature R and the five-form 
square F| will suffice. Note that the five-form is always accompanied by a power of 
as is the case for RR forms. 
The function F in fl3.1|) will be taken arbitrary at this stage. Notice also that we 
did not add non-linear terms in the kinetic term of the dilaton, in the string frame. 
Following many works and explicit analysis of the string cr-model up to four-loops it 
was conjectured that there is only a linear term in {dip) 2 . We will assume this 
to be true here. 



3 In particular the holographic conformal anomaly |3(| and the shear viscosity |37j can depart 
from their universal values only if higher derivative structures involving the Riemann tensor are 
present. 
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For weak curvatures we can expand 

F = ~ + R+^ + 0(R 2 ,R£,e) , 5c=10-5 = 5 (3.2) 

to obtain the standard two derivative tree effective action including the dilaton po- 
tential (the first term) present due to the fact that we are in non-critical string 
theory. 

As the four-form potential is non-propagating in five dimensions, it is appropriate 
to "integrate it out". To do this we derive its equations of motion, solve them, we 
then substitute back into the equations of motion of the others fields (metric and 
dilaton) and we find the new action from which these equations stem by variation. 

Varying (Of) with respect to the four-form we obtain the equation 



V (F ? F^) = (3.3) 
where F$ = 9F ^^> , The solution of (|3.3|) is 

N e 

TP TP — c * i VP (TT (n A\ 

r € 1 tivpoT — „ 2 /t: \°-^) 

c Ads yy 

where called the constant N c and inserted the AdS scale ^Ads that will be introduced 
later in ( |4.9| ) to make N c dimensionless 4 . All we know of course is that this constant 
is linear in the number of colors, but unlike the critical M = 4 case we do not know 
the precise coefficient. To obtain this we need to know the D% solutions in the non- 
critical theory in question and their tensions. We will still call the flux constant iV c 
though from now on, keeping in mind that this is only proportional to the number of 
colors. What is important is that we are working in the limit where this flux is sent 
to infinity keeping 

A = N c e* (3.5) 

fixed. Indeed this combination is proportional to the 't Hooft coupling of YM, at 
least in the UV. Squaring ( |3.4| ) we obtain 

e AdS 

This is an algebraic equation that involves, £, R, A, and we must solve it implicitly to 
obtain £ as a function of R, A. This can be used to obtain the following differential 
formulae, 

(F e + 2£F se )d£ + 2£F^di?=2£F^ , j| = - ^ (3.7) 



4 Although the AdS length will emerge later, it is inserted here for economy, in order to make 
later formulae simpler. 
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useful for the variations of the action with respect to the other fields. The equations 
for the dilaton and graviton read 

n<f>-{d<f>) 2 + I(F-£F e ) = (3.8) 

1 e 2 * F 2 

F R + e 2<t> (g^U - V M V,)(e- 2< ^) + Ad^d u <P - -(4(«90) 2 + F)g flv - i^^p 1 = 

(3-9) 

where Fr = 9F qr^ ■ Substituting from ( |3.6[ ) into ( |3.9| ) we obtain the equation 



Fr + e^ig^D - V^V^e" 2 ^) + 4«9 M <R0 - ^(4(<90) 2 + F)g, w + f F^ u = 

(3.10) 

Equations ( |3.8| ) and ( |3.10| ) can be obtained from the following equivalent action 

Stree = / (fx^e^ [4(«90) 2 + F(R, £) - 2^(R, 0] , (3.11) 

where in this action £ = £,(R, A) is an algebraic solution of ( |3.6| ). Note that this 
action is 0(N 2 ) if we use variables that are finite in the large- N c limit as follows 



S tme = M 3 N 2 J d'x^g- 



' f)\2 

4— + F(R,t)-2tFs(R,0 



(3.12) 



The conclusion of this analysis is that in the end of the day we must solve the 
following two equations subject to the algebraic condition ( |3.6| ) 



Fr R, u + e^ig^D - V,V u )(e~ 2 *F R ) + 4^<R0 - ^(4(<90) 2 + F)g fiu + £ = 

(3.13) 

4D</> - 4(d0) 2 + F - ^Ff. = (3.14) 

4. The UV regime 

It was already argued in section that the most reasonable description of the asymp- 
totic UV geometry is as an AdSs near-boundary region. This description geometrizes 
the asymptotic conformal invariance, and brings us in line with what was understood 
in the best studied case of Af = 4 superYM . 

The conformal invariance in QCD however comes together with a coupling con- 
stant that vanishes in the UV. Therefore we expect that as we approach the boundary 
A — > 0. In the metric ansatz ( |2.23|) our expectation translates into the fact that near 



the boundary 

b(r) - ^ (4.1) 
r 

the AdSs warp-factor in Poincare coordinates. Moreover in this regime the intuition 
is similar to the M = 4 case: r is serving as the inverse of the energy scale. 
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QCD tells us that near the UV, the 't Hooft coupling depends on the energy as 

^ = 6 log| + O(loglog|) (4.2) 

Using r as the inverse of the energy we deduce that the solution for the 't Hooft 
coupling A near the AdS boundary must look like 

— = -b log rA + C(log log rA) (4.3) 
A 

and therefore A — > as lo ^ rA as we approach the boundary r — > 0. 

We conclude that asymptotically close to the boundary R — > R* = — -pr- and 
A -> 0. 

Consider now equation (|3.6|) . Since the right-hand side vanishes near the bound- 
ary, the left hand side must vanish also. This can happen in two possible ways: as 
A — > either £ — > or £ — > ^ 0, with i^(£*) = 0. We will examine both options 
in turn and show that only the second one can consistently happen. 

4.1 Vanishing £ in the UV 

If £ — > 0, as A — > then £ must vanish as 5 

F^(R*, 0) £^ dS 

Assuming an AdS 5 (constant curvature) solution R = equations (|3.11|) - (|3.l4|) 
imply to leading order (A = 0) that 

F(R m ,0) = F R (R*,0) = (4.5) 

This in turn implies that the AdS curvature R* must be a double root of F(R, 0). 

We now move to the next order and perturb the leading solution R = i?*, A = 
to R = R* + 5R and non-zero but small A, 



A 2 



F^(R*, 0) £\ dS 



Fk(r.,o) a 2 ,,/■•,/,■(/,'.■()) 

F|(^,0)£i ds F € (i2.,0) + ' 



(4.6) 



A 2 A 2 
F ^F e (^,0) ^ dS + "' ' ^^F c (^,0)£i ds + "' (47) 

where we have used (|4.5|). To investigate this next order solution we use the desired 
asymptotic form of the dilaton 

A ~ — : - —r-. . + • ■ ■ — ► <f) ~ constant — log[— 6 log(Ar)] + ■ • ■ (4.8) 
6 log(Ar) 



5 This assumes that remains constant as £ — > 0. Otherwise it can be shown that there is no 
solution. 
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and using the AdS metric 



P 2 2fl 
ds 2 = ^f-[dr l + dx^dx ll \ , R* = - r (4.9) 

r M.dS 



to leading order we obtain 



4 



□0 ~ n , - — + • • • , (d0) 2 ~ tf2 : .. , 9 + • • • (4.10) 
£ 2 dS log(Ar) YJ £i dS log(Ar) 2 

Equations (|4.7|) and ( |4.9j ) are now incompatible with ( |3.14| ). Therefore, the 
starting assumption cannot be correct for the string theory we are seeking 6 . 

4.1.1 Non- vanishing £ in the UV 

In this case the function F% must have a zero as a function of £ at a non-trivial value 
We can therefore parametrize for convenience 

F ~ co(i?) + ^(£ - ^(i?)) 2 + O [(£ - ^(i?)) 3 ] (4.11) 



and obtain 7 

£ = + 5£ ~ ± — A == + 0(X 2 ) (4.12) 

Ci(-R) «AdS v£*(-R) 

We keep both signs here but below we will see that only the minus sign is relevant. 

Again the gravitational equation implies that for AdS to be the leading solution 
(at A = 0) we must have 

c °( R > )=0 • mL,=° < 4 - 13) 

Using the above we obtain that F is zero to next order and the first non-trivial 
contribution is at quadratic order 



1 a 2 Cp , 

2ciCR0 £ 2 AdS UR*) + 2dR 2 



F (R, = o ^2 ,, n ^ + n^(R*)(R - R*) 2 + ■■■ (4-14) 



and therefore subleading while 



A 



{•AdS 



6 A more exotic possibility is to cancel the offending term against an SR 2 term from the next 
order expansion of F. This however will give an effective expansion in powers of A^ and not integral 
powers of A, a fact at odds with QCD perturbation theory. 

7 For similar reasons the more general possibility F ~ Cq(R) + — £*(-R))° +1 + ■ ■ ■ with 

a > would imply A ~ log(Ar) - ". Therefore only the case a = 1 is relevant for YM. 



- 19 - 



Now it is possible to solve ( 3.14 ) to leading order and match 



^AdS 



\l (R* 



16 



(4.16) 



To have asymptotic freedom we must choose the minus sign in (|4.12| ). 

Continuing further to the trace of the gravitational equation ( p,13|) we obtain 



F R R + 4DF R - 16d(j)dF R - 8F R n<j> + 10(<9</>) 2 - -F + 5£F C = 



(4.17) 



Taking into account that to leading order F R = 0, to next order F = and ( |4.13j ) 
we obtain that ( [4.171) becomes to next order, 



(4D + R*)5R 



where we have used the fact DA 



5 + 



4 

" b £l dS log(Ar)2 



(4.18) 



-AdS 



which gives a subleading contribu- 



tion. This equation gives the following leading modification to the AdSs metric 



b= - 

r 



1 + 



w 



with 



log(Ar) 



w 



+ 



5R 



AOw 



log(Ar' 



+ 



- 5 + t^ R *UK 



(4.19) 



(4.20) 



cg(i?,) 80 R, 

Let us pause and review what we have found. Near the boundary the leading 
solution is AdSs and the dilaton is such that A = 0. 



As all dilaton-dependent parts of the effective action ( |3.11| ) are subleading 
near the boundary, the leading AdSs solution must be supported by curvature 
alone. This is the essence of the conditions ( |4.13j ). In particular the non-critical 
dilaton potential and the corrections coming from integrating out the four-form 
are subleading near the boundary due to asymptotic freedom. 

The fact that the AdS solution is supported by curvature alone implies that 
modulo accidents, the AdS curvature scale ^Ads is of the same order of mag- 
nitude as the the string scale £ s . Do we know string backgrounds with such 
a property? In a sense yes, although the backgrounds we know are somewhat 



simpler. In particular coset CFTs, |39j share some features with what should 
happen here. In coset CFTs there is a special frame in which the solution to 
the a model conditions can be thought as exact, but in other schemes it obtains 
corrections. However, unlike what is expected to happen here, in coset models 
one can vary (in a discrete fashion) the curvatures by varying the levels of the 
current algebra. In a sense we must have a solution where the curvature cannot 
be varied. 
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• Unlike situations in critical string theory, the asymptotic AdS geometry here is 
not supported by RR flux. Therefore, to leading order near the boundary the 
a-model is "conventional. It is to next order that the flux enters the solution 
and the coupling starts to run. 

• So far we have seen the importance of non-linearity of the curvature part of the 
effective string action in order to find an leading AdS solution. Note however 
that in order for the function F% to have a root at a non-zero value of £, the 
function F(£) must also be non-linear. Therefore, we also need the higher 
derivative corrections of the four- form to find the asymptotic solution of YM. 

• Many authors have suggested that the key to understanding the UV region 
of the QCD string theory is a theory of an infinite number of massless higher 
spin fields. The idea behind is coming from M = 4 holographic intuition. 
There, A — > implies that a' — > oo and the whole string spectrum is becoming 
massless. We see here that things work differently. The string has finite stiffness 
at the boundary, although it is soft, with a fundamental string tension that is of 
the same order of magnitude as the background curvature. On the other hand, 
all three and higher-point connected correlators will vanish near the boundary 
for a simple reason. When properly normalized, they are all multiplied with 
positive powers of A that vanishes near the boundary. 

• In this section we worked out explicitly the first non-trivial order of the string 
equations. We have "imposed" both an asymptotic AdS solution and a leading 
running of the coupling constant. For the rest, assuming genericity, there is 
a regular perturbation theory in inverse powers of logs, that is similar to the 
perturbative expansion in perturbative QCD. The various coefficients arise from 
the expansion coefficients of the non-linear string effective action around the 
vacuum solution. Without complete control of the non-critical string theory 
they cannot be calculated. Despite this, the structure of the near-boundary 
perturbation theory is clear. Moreover the simplifying assumptions we made 
about the effective action do not seem to modify the conclusions above. 

5. The IR regime 

We have seen that the general structure of the string effective action and some 
simple assumptions on the UV asymptotics indicate the presence of the standard 
YM perturbation theory in the UV. The situation in the IR is much murkier for a 
the simple reason, that no guiding principle like perturbation theory is known to 
exist. We do have specific expectations however from the strong coupling region 
of QCD. In particular we expect confinement and a discrete and gapped glueball 
spectrum. However these requirements are fairly indirect to guide us in the IR. 
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A first question we would like to ponder is what happens to the dilaton in the 
IR. The most natural expectation is that increases without bound, so that A — > oo. 
There have been minority claims of an IR fixed point in QCD but such claims are not 
in our opinion credible. There is also the alternative that the coupling asymptotes 
to a finite large value in the IR. Although this is not excluded, we will not entertain 
it here for two reasons. The first is that we did not find a good way to implement 
this possibility while keeping all the properties we expect from YM in the IR, in 
particular confinement. The second is that even if A — > oo in the IR for most 
observables , there may be a maximum finite value for the 't Hooft coupling if all 
low lying wave-functions have support away from the A — > oo limit. 

Vacua with a runaway dilaton are known in string theory. The simplest is the 
linear dilaton vacuum that is simple enough to define, but has consistently puzzled 
researchers for two decades. Due to advances in the understanding of Liouville theory, 
the associated matrix models and more recent advances in the associated world-sheet 
CFTs we have today a fairly good idea of the physics in such a background. A 
crucial ingredient is that there should be a sufficient screening of the strong coupling 
singularity, that in Liouville theory is achieved via the "Liouville wall". We will 
see later on that a linear dilaton background in the IR is marginally compatible 
with what we expect from YM at strong coupling. In particular we can show that 
we have confinement and a mass gap, but the spectrum is continuous above the 
gap. However, by slightly modifying the background we will obtain also a discrete 
spectrum as well as linear asymptotic trajectories. Moreover as will see, there is 
a sense in which the strong coupling singularity is screened: all local low-energy 
observables do not get contributions from arbitrarily close to the singularity. This is 
what we will call a "repulsive" singularity that is a more constrained concept than 
the "good" singularities of Gubser, ||23|| . 

Another intuition is emerging from the M = 4 paradigm of AdS/CFT: at strong 
coupling we could expect a good effective description in terms of a two-derivative 
action. Although this is rather transparent in standard AdS/CFT, it is less clear 
here, because the 't Hooft coupling is not constant. We will however take it as a 
principle and we will see how far we can go. An important ingredient in order to 
implement this idea is that the curvature in the string frame must be very small in 
the IR. This ties well together with the linear dilaton paradigm as in that case the 
curvature in the string frame vanishes. We would therefore investigate the possibility 
that there is a vacuum solution in the IR with A — > oo, a small curvature in the string 
frame, and subleading contributions from higher derivative terms. 

To investigate this we will expand now the string effective action in ( |3.12| ) in 
powers of the curvature, remembering that £ is an implicit function of R and A from 
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F(R, - 2^(R, = E ^n(A) i? n 
so that the action becomes 



n=0 



Stree = J d*Xy/ge 



2,p 



4(90) 2 + z (\) + z 1 (x)r + J2 z n {>) R n 



n=2 



(5.1) 



(5.2) 



We now make a conformal transformation 

g»u^f(4>)g,v , R^j[R-4n\ogf-3(d\ogf) 2 } , D^j 



with 

to obtain 



Stree — / d Xy/g 



R 



4-3 



□ + -V^(log/)V M 
(5.3) 

(5.4) 

ef Zn A Z n (A)Zf (2 - 5) 

— — ^ + ) 



w) 2 + ^ + E 

n=2 



We may now define a new scalar $ with a canonical kinetic term 



d$ / 4 3 (1+^) 



(5.5) 



(5.6) 



Asking for positivity inside the square root, severely constraints on how Z\ depends 
on (p. In particular we find that 



Zi((j>) = C\e^ + subleading , as <f> — > 00. 
The action (|5.5|) can now be written as 

4 



(5.7) 



Stree = I d 5 X^/g 



R- (d$) 2 + V ($) + ]TK($) (R + y 



n=2 



with 



i(2n-5) 



e |*(»-i) 



(5- 



(5.9) 



4 T 

We may now return to a "string frame" for $ as g^ u — > e~z to obtain 



5ms = / d b x^/ge 



-2$ 



i2 + 4(«9$) 2 + W ($) + nn 



n=2 



(5.10) 
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where 



Wb(*) = |?4r , = ^ F 2J ^ , F = Zl e^ (5.11) 

TZ = R+^DlogF -UdlogF) 2 (5.12) 

From now on we will always call A = e* the 't Hooft coupling and use it interchange- 
ably with $. 

Before we proceed further we should discuss our expectations on the dependence 
of the coefficients Z n on the 't Hooft coupling in the IR. We do expect the leading 
dependence to be the same for all Z n . Moreover, we expect that their growth in 
the IR is bounded. This is already suggested by the positivity bound on Z\ that we 
pointed out in ( |5.7|) . 

In Jl8| an analysis of various dilaton potentials was performed, assuming a two- 



derivative action, in order to find which ones give properties that we expect from 
YM, namely, confinement, mass gap and discrete spectrum. At a two-derivative level, 
positive energy implies that all Lorentz invariant ansatze lead to an IR singularity 
(the only other alternative is an IR AdS space). Therefore part of the criteria we 
used included that this singularity is repulsive 8 : spectra can be computed without 
the need of extra boundary conditions at the singularity and string world-sheets 
should not come very close to the singularity. 

It was found that potentials fall in several categories, that we parametrize their 
asymptotics as V ~ A^(logA) p . The presence of the log terms matter only when 



Q - 4 



3 • 

Potentials with Q < | do not confine. 

Potentials of the form V ~ A 4 / 3 (log A) 2 ^ 1 , a > 0, satisfy all criteria, and 
moreover they have the property that the string frame curvature vanishes in 
the IR. For a = 2 they generate asymptotically linear trajectories, a = 1 is the 
marginal case of the linear dilaton vacuum. This case does not have a discrete 
spectrum. 

Potentials with Q = 4/3 and P > 1, confine properly, but the string frame 
curvature blows up at the singularity 

Potentials with | < Q < 4 -^ confine properly but the string frame curvature 
blows up at the singularity. 

Potentials with Q > have a "bad" (non-repulsive) IR singularity. 



^This is a stronger condition than the one in Gubser's classification, 
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In |18| a detailed analysis of the glueball spectra of different confining potentials 
were performed. It was found that only the "soft potentials" V ~ A 4 / 3 (log A) , a > 
are capable of giving spectra that are reasonably close to the lattice glueball spectra 
All other potentials produced glueball splittings (after adjusting/fitting parameters) 
that cannot be accommodated by the lattice data. 

In view of this we may analyze the possibilities in the IR using the actions ( |5.S| ) 
and ( 5 . 1 (J| ) and the principles that Z n jZ\ — ■> constant as A — > oo and that the bound 



applies. 

We find the following possibilities 



4 , 



1. When Z\ ~ we find a potential that leads to non-confining behavior, and 
higher derivative corrections that are no controlable. 

2. When Z\ ~ e b ^(j) c with b < | leads to non-controlable higher derivatives with 
the exception of b = —4 + 2^/7 ~ 1.29. In the latter case the potential is 
Vq ~ es $3 C and this is an good confining potential with vanishing string 
frame curvature provided c < |. In this case all higher derivative corrections 
are subleading. 

3. When Z\ ~ constant plus subleading we end up with a non-controlable higher 
derivative behavior with the exception of Z\ — 1 — |r + • • • . In that case the 

leading effective potential is V ~ e^®Q c which gives rise to confining behav- 
ior and vanishing string frame curvature. In this case the higher derivative 
corrections are again subleading. 

4. Any other asymptotic behavior leads to non-controlable higher- derivative cor- 
rections. 

It is not obvious which of the two favorable possibilities found above are realized 
in QCD. We cannot even exclude the cases that have non-controlable higher deriva- 
tive corrections. It is noteworthy though that both favorable cases 2 and 3, give the 
same order of magnitude corrections to the IR running of A. Starting from the Ein- 
stein frame effective action and keeping the leading quadratic term (that is □ log F) 
we can directly estimate that its influence on the leading order solution scales as 

A (log A) 2 



4 



for all soft potentials V ~ AaA - ^ - , a > 1. 

With a bit of optimism and guided by the intuition above we will accept that the 
IR of QCD is governed by a soft potential, that once taking the existence of linear 
trajectories into account, it should have the following asymptotics V ~ X^y/\og\. 
Moreover, in such the higher-derivative corrections are suppressed. 
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5.1 The 't Hooft coupling in the IR 



We have already seen that in the IR, the canonically normalized dilaton field $ is 
non-linearly related to the original string theory dilaton <fi. In particular in case 2 of 
this section with b = —4 + 2\/7 this relation is 

$ = (3- Vl)<p + c\og<p + ■ ■ ■ (5.14) 

where c a constant, while for case 3 

$ = + clog0H (5.15) 

where c is again a constant. 

There is however another source of corrections to the identification of the 't Hooft 
coupling, [ET7| if it defined using the D% brane world- volume action 9 . The general 



form of the kinetic term for the gauge fields on a D 3 brane is expected to be 

C F 2 = e~*Z(R, i) Tr[F 2 ] (5.16) 

where Z(R,£) is an (unknown) function of curvatures and the five-form field strength. 
At weak background fields Z ~ — | + • • • . 

In the UV regime, using the formulae of section [4.1.1| we obtain 



C F 2 = N c Tr[F 2 ] - 



Z (i?*, £*) 



A 



from which we can identify the QCD 't Hooft coupling as 

A ^(i?*, £*) A 2 



+ 0(A 2 



(5.17) 



A 



QCD 



+ 0(X 2 



(5.U 



indicating that there are "perturbative" corrections to the identification to the 't 
Hooft coupling. In the IR, the relation is more complicated and to leading order in 
the string curvature it becomes 



A 



A 



QCD 



Z(0,U0, A) 



(5.19) 



We do not know what this non-trivial function of the dilaton is. We should stress 
though that this is a particular definition (scheme) for the 't Hooft coupling constant 
(that we call the D 3 scheme, and it may be very different (especially in the IR) from 



other schemes used in lattice calculations, 46 



3 We thank K. Kajantie whose question prompted this discussion. 
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6. Schemes and scheme-dependence 

As mentioned in the previous section, there are several choices in the theory that 
amount to a choice of scheme. Schemes in QFT correspond to coordinate choices 
in the space of couplings. One can parametrize the space of couplings in different 
fashions but that is not affecting physical observables that should "coordinate invari- 
ant". A particular scheme choice is that of the renormalization group scale used in 
renormalizing the theory. Here again physical observables should be RG-invariant. 

In the holographic context scheme dependence related to coupling redefinitions 
translates into field redefinitions for the bulk fields. As the bulk theory is on-shell, 
all on-shell observables (that are evaluated at the single boundary of spacetime) are 
independent of the field redefinitions showing the scheme-independence is expected. 
Invariance under radial reparametrizations of scalar bulk invariants is equivalent 
to RG invariance. Because of renormalization effects, the boundary is typically 
shifted and in this case field redefinitions must be combined with appropriate radial 
diffeomorphisms that amount to RG-transformations. 

In particular, the holographic definition of the "Energy" should not affect physi- 
cal quantities, provided it is reasonable: (a) it is monotonic and (b) it vanishes at the 
ultimate IR. Of course quantities like the /3-function do depend on the definition of 
the energy (as well as the definition of the coupling). This is expected as a /3-function 
is not an invariant (physical) but a vector in coupling constant space, as well as in 
energy space and therefore changes in different schemes. 

7. A phenomenological model 

In view of the discussion in the previous sections, we should ask the question: to 
what extend a tractable simplified model can be devised that is reasonably close to 
what we expect to happen in the low-energy theory. 

The answer to this question in the IR is straightforward. We have argued that 
the leading IR asymptotics are described by a dilaton gravity system with a soft 
potential, and no higher-derivative corrections. 10 

Vo ~ A a \/log\ as A — ► oo 

(7.1) 

The non-trivial issue is the UV. There we have seen that higher-derivative cor- 
rections, especially in the curvature, are essential for the appearance of the asymp- 
totically AdSs solution and a logarithmically- vanishing coupling constant. Imple- 
menting this at the present stage seems out of reach if we insist in a model that we 



S 



HQCD 



R 



10 Dilaton gravity with a potential was also argued to have similar thermodynamic properties to 
QCD in g§. 



can calculate observables. There is short-cut though that has the advantage that it 
can capture some of the salient features of the UV solution and dynamics while it 
remains manageable. As expected, it has also shortcomings that we will list below. 

The shortcut is to assume a similar action as in ( [7.1]) but with a potential that 
has different weak coupling asymptotics. 

Vo(X) = ^- [1 + Vl X + v 2 X 2 + v 3 X 3 + ■ ■ ■} as A^O (7.2) 

^AdS 

This form is indicated by the fact that if v% (that implies that the ultimate 
boundary is not a minimum of the potential) and the dilaton boundary conditions are 
carefully chosen (as detailed in |L8| pOR), the vacuum solution for the 't Hooft coupling 
and the metric has the same structure in terms of the expansion in inverse logs, as 
the one advocated here in section f| and which agrees with QCD perturbation theory. 
In particular once we chose a definition of the energy scale 11 , then the dimensionless 
coefficients v n are in one to one correspondence with the /3-function coefficients. 

Moreover the basic property that as we reach the boundary three- and higher- 
point connected correlators vanish is satisfied, since A — > 0. 

There are further consistency checks that such an approach in the UV captures 
quite a few important properties of the system. We list them below 

• It captures the correct number of UV boundary conditions in accordance with 
YM 

• It gives the correct perturbative running of the coupling once the analogue 
of Aqcd is set. This was verified in [JTSJ by reproduced the experimentally 
measured value of the strong coupling constant at E — 1.2 GeV, using as an 
input the lowest glueball mass. 

• It provides the correct UV boundary conditions for the physical fluctuations of 
the system, the glueballs, p~8f| . 

• It is compatible with the existence of a finite temperature phase with the correct 



physical properties, [19, 20 



• It provides the correct free-field limit at high temperature and its first loga- 
rithmic correction. 

There are also issues that are not as they should be. They include 



11 In JTt], |l^| it was argued that a good (but not unique) choice is the scale factor of the metric in 
the Einstein frame. Moreover different choices that are perturbatively related affect the /3-function 
coefficients beyond the two-loop level. 
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The conformal anomaly. It is well known [3(J that any theory with a leading 



two-derivative bulk description must be such that the two anomaly coefficients 
are equal a = c to leading order, 0(N%). On the other hand, QCD with or 
without quarks, has a ^ c to leading order in N c 



12 

' c 



The shear viscosity. It was shown in [40[ that two-derivative theories of gravity 



coupled to matter, have all the minimum constant value of the 2 ratio. For the 
IR regime, as we argued earlier this is probably a very good approximation. 
However as we approach the UV, this fact is more and more at odds with what 
happens in perturbative QCD where - diverges as the coupling vanishes. 

In the UV there are some spurious extra logs in some quantities. They reflect 
the fact that in the string theory dual of YM, the asymptotic AdS metric is 
achieved in the string frame and not the Einstein frame. This is a non-trivial 
issue and can be ascertained from various points of view. The analysis in 
section |3] is performed in the string frame. It would not go through as such if 
an AdS metric in the Einstein frame is advocated. On the other hand, for the 
phenomenological model discussed in this section, it is not consistent with the 
equations of motion to impose the AdS metric in the string frame. We therefore 
impose it on the Einstein frame. The difference is an extra multiplicative factor 
of (log(Ar))3 since the dilaton is asymptotically logarithmic. The effects of 
this result in mild differences in some quantities. For example a calculation 
of the short distance asymptotics of the static quark potential produces an 
an extra log factor of the distance, multiplying the 1/r of the Coulomb law. 
This modification is mild and was shown to fit as well baryonium data as the 
standard Cornell potential, p£2fl . 



• Higher order correlators of the basic operators of the theory, or higher per- 
turbative (ie. log corrections) to low-order correlators may be quantitatively 
different although generically they remain qualitatively correct. 

7.1 Parameters, fits and predictions 

The action (|7.1] j contains as parameters two scales, M p and ^Aas from the potential 
( |7.2p as well as in principle many dimensionless parameters hidden in the potential. 
There is also the string scale £ s that enters the fundamental string action. We do 
not know ab initio how such scales are related. However we may observe that Iaos 
is not really a physical parameter by a choice of energy scale. 

Moreover there is dependence on boundary conditions for the vacuum solution. 
There are three boundary conditions needed for the gravity-dilaton system. One of 



12 It can be shown the for CFTs, with a minimum of Af = 1 supersymmetry only fundamental 
matter can contribute to a — c, JL3| . For theories with a weak coupling limit, conformal invariance 
implies that only theories with the Af = 4 spectrum have a — c. 
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them is a gauge artifact, as it can be identified with a translation of origin in the radial 
coordinate. Another is fixed by demanding that correct perturbative asymptotics of 
the dilaton and the absence of a bad singularity in the bulk. The third corresponds 
to Aqcd- If is defined in a reparametrization invariant way as 



A 



_1 lim 



6(A)- 





1 


exp 


b \_ 



(7.3) 



where 6(A) is the scale factor of the Einstein metric in (|2.23j ) written as a function of 
A, and 6 , 61 are the one and two-loop coefficients of the perturbative YM /3-function. 
A is fixed once we specify the value of the scale factor 6(A) at a given Ao- 

M p does not affect spectra of particles but it affects interactions. It also affects 
the size of the free-energy at finite temperature. It can be fixed, 



{M p £ 



AdSy 



45 7T Z 



(7.4) 



by demanding that the free energy for T — > oo asymptotes to that of a free gas of 
gluons. The fundamental string scale can be fixed by calculating the inter-quark 
potential and comparing with the associated lattice results for the effective string 
tension. 

In practice the rest of hidden parameters of the potential are truncated to a small 



set. In |[L8| a one-parameter fit could reproduce all known lattice glueball spectra with 
values inside the error bars, and all the characteristic features of the thermodynamics, 
OJ. Adding a second parameter the full equilibrium thermodynamics fits perfectly 



well, 44 



Once the parameters have been fixed, then one can calculate other quantities. 
They include the full spectrum two and higher-point functions, the full equilib- 
rium thermodynamic quantities ||44|| , as well as transport data. 
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